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In this paper we present AfRne.m - program for computations in represen- 
tation theory of finite-dimensional and affine Lie algebras and describe imple- 
mented algorithms. Algorithms are based upon the properties of weights and 
Weyl symmetry. The most important problems for us are the ones, concern- 
ing computation of weight multiplicities in irreducible and Verma modules, 
branching of representations and tensor product decomposition. These prob- 
lems have numerous applications in physics and we provide some examples 
of these applications. The program is implemented in popular computer al- 

■^ gebra system Mathematica and works with finite-dimensional and affine Lie 

t^^ algebras. 
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Nature of problem: 

Representation theory of finite-dimensional Lie algebras has many applications 
in different branches of physics, including elementary particle physics, molecular 
physics, nuclear physics. Representations of affine Lie algebras appear in string 
theories and two-dimensional conformal field theory used for the description of 
critical phenomena in two-dimensional systems. Also Lie symmetries play major 
role in study of quantum integrable systems. 
Solution method: 

We work with weights and roots of finite-dimensional and affine Lie algebras and 
use Weyl symmetry extensively. Central problems which are the computations of 
weight multiplicities, branching and fusion coefficients are solved using one general 
recurrent algorithm based on generalization of Weyl character formula. We also 
offer alternative implementation based on Freudenthal multiplicity formula which 
can be faster in some cases. 
Restrictions: 

Computational complexity grows fast with the rank of algebra, so computations 
for algebra of rank greater than 8 are non practical. 
Unusual features: 

We offer the possibility to use traditional mathematical notation for objects in 
representation theory of Lie algebras in computations if AfRne.m is used in Math- 
ematica notebook interface. 
Running time: 

From seconds to days depending on rank of algebra and complexity of representa- 
tions. 



1. Introduction 

Representation theory of Lie algebras is of central importance for different 
areas of physics and mathematics. In physics Lie algebras are used for the 



description of symmetries of quantum and classical systems. Computational 
methods in representation theory have a long history pp , there exist numerous 
software packages for computations related to Lie algebras [2], [3], [HE], [6], 

m- 

Most popular programs [2], [3], [H], [S] are created to study representation 
theory of simple finite-dimensional Lie algebras. The main computational 
problems are the following: 

1. Construction of root system which determine the properties of Lie al- 
gebra including its commutation relations. 

2. Weyl group traversal which is important due to Weyl symmetry of root 
system and characters of representations. 

3. Calculation of weight multiplicities, branching and fusion coefficients, 
which are essential for construction and study of representations. 

There are well-known algorithms for these tasks [8], |9], [1], [10]. The third 
problem is most computation-intensive. There are two different recurrent 
algorithms which are based on Weyl character formula and Freudenthal mul- 
tiplicity formula. In this paper we analyze them. 

Infinite-dimensional Lie algebras also have growing number of applica- 
tions in physics for example in conformal field theory and study of integrable 
systems. But infinite-dimensional algebras are much harder to study and 
number of available computer programs is much smaller. 

Affine Lie algebras [TT] constitute important and tractable class of infinite- 
dimensional Lie algebras. They are constructed as the central extensions of 
loop algebras of (semi-simple) finite-dimensional Lie algebras and appear nat- 
urally in the study of Wess-Zumino-Witten and coset models of conformal 
field theory [H], [13], [H], [IS]- 

The structure of affine Lie algebras allows to extend computational al- 
gorithms created for finite-dimensional Lie algebras [7], [16], [17]. The book 
[T7] with the tables of multiplicities and other computed characteristics of 
affine Lie algebras and representations was published in 1990. But we are not 
aware of software packages for popular computer algebra systems which can 
be used to extend these results. We address this issue and present AfRne.m 
- Mathematica package for computations in representation theory of affine 
and finite-dimensional Lie algebras. We describe the features and limitations 
of the package in present paper. We also provide representation-theoretical 
background of implemented algorithms and present examples of computa- 
tions relevant to physics. 



The paper starts with an overview of Lie algebras and their representation 
theory (Sec. [2]). Then we describe datastructures of AfRne.m used to present 
different objects related to Lie algebras and representations (Sec. |3]) and 
discuss implemented algorithms (Sec. 111). Next section consists of physically 
interesting examples (Sec. [s]). The paper is concluded with the discussion of 
possible extensions and refinements (Sec. |6]). 

2. Theoretical background 

In this section we remind necessary definitions and present formulae used 
in computations. 

2.1. Lie algebras of finite and affine types 

Lie algebra g is a vector space with bilinear operation [-, ■] : ® — )■ 0, 
which is called Lie bracket. If we choose the basis Xj in q we can specify 
commutation relations by the structure constants Cijk- 

[X\X^]="^&^X^ (1) 

k 

Lie algebra is simple if it contains no non-trivial ideals with respect to com- 
mutator. Semisimple Lie algebra is a direct sum of simple Lie algebras. In 
present paper we treat simple and semisimple Lie algebras. 

Maximal commutative subalgebra {Cartan subalgebra) of g is denoted by 
[)g. We denote the elements of basis of f)g by H\ 

Killing form on q gives a non-degenerate bilinear form (-, ■) on f)g which 
can be used to identify f)g with the subspace of the dual space [)* of linear 
functionals on [)g. Weights are the elements of [)* and are denoted by Greek 
letters /i, i/, cu, A . . . 

Special choice of basis gives compact description of commutation relations 
(Q. This basis can be encoded by the root system which is discussed in 



Section 2.2 (See also [I 



Loop algebra Lg = g ® C[t,t ^], corresponding to semisimple Lie algebra 
Q, has commutation relations 

[X*r, XH""] = t"+'" J2 '^k^^ (2) 



Central extension leads to the appearance of additional term 

[X'r + ac, X^'r + Pc] = r+'" Yl '^k^' + (^'' ^'>^n+m,oc (3) 

k 

This algebra g = g® C[t,t^^] © Cc is called non-twisted affine Lie algebra 

[m, [SQlEI], [I7J. 

We do not treat twisted affine Lie algebras in present paper. 

2.2. Modules, weights and roots 

Let Q be finite-dimensional or affine Lie algebra. 

Then g-module is a vector space V together with a bilinear map qxV ^ 
V such that 

[x,y] ■ V = X ■ {y ■ v) — y ■ {x ■ v), ior x,y G q,v E V (4) 

Representation of algebra on a vector space V is the homomorphism g — >■ 
gl{V) from to Lie algebra of endomorphisms on vector space V with the 
commutator as the bracket. 

For an arbitrary representation it is possible to diagonalize the operators 
corresponding to Cartan generators W simultaneously by special choice of 
basis {vj} in V: 

W . V, = v]v, (5) 

Eigenvalues i/* of Cartan generators on element of basis Vj determine weight 
Uj G f)* such that i'j{H^) = i^j- Vector f G V^ is called weight vector of 
weight A if Hv = \j{H)v, WH G f). Weight subspace consists of all weight 
vectors Vx = {v e V : Hv = Xj{H)v, \/H G f)}. Weight multiphcity 
m\ = mult (A) = dimVA is the dimension of weight subspace. 

The structure of module is determined by the set of weights since the 
action of generators E'^ on a weight vectors is 

E"" -VxOC Vx+a (6) 

Module structure can be encoded by the formal character of module 

chV" = ^ nixc^ (7) 

A 

Character chV" G £^ is an element of algebra S generated by formal exponents 
of weights. Character can be specialized by taking its value on some element 
eoff). 



Lie algebra is its own module. The action of generators is called adjoint 
and it is given by the bracket adxY = [X, Y]. Roots are weights of adjoint 
representation of q. They encode the commutation relations of algebra in 
the following way. Denote by A the set of roots. For each a G A there exist 
a root —a G A and generators E'^, E~'^ such that 

[i7\E"]=a*E" (8) 

if a + /3 G A 




[^"'^T = <!raE.«^^\ if« = -/3 (9) 

otherwise 

Given a root system A we can choose a set of positive roots. This is a 
subset A"*" C A such that for each root a G A exactly one of the roots a, —a 
is contained in A"*" and for any two distinct positive roots a,/3 G A+ such 
that a + /3 G A their sum is also positive a + /3 G A"*". Elements of — A+ are 
called negative roots. 

A positive root is simple if it cannot be written as the sum of positive 
roots. The set of simple roots $ = {o-i} is a basis of f)g and each root 
can be written as a = ^^ niUi with all nj non-negative or non-positive. In 
case of finite-dimensional Lie algebra g simple roots are numbered from 1 to 
rank of the algebra i = 1, . . . , r, r = rank(0). By numbering simple roots 
with an index i we introduce lexicographic ordering on the root system A. 
Highest root with respect to this ordering is denoted hy 6 = X]j=i r'^i^ii 



coefficients Oj are called marks. It is also highest weight (See section 2.3) of 
adjoint module. Comarcs are equal to a^ = ^"'^"'-^ aj. 

Although the full set of roots A is infinite for affine Lie algebra q the set 
of simple roots $ is finite and its elements are denoted by Oq, . . .Or where 
r = rank(g). The roots ai, ...,«,. are the roots of the underlying finite- 

o 

dimensional Lie algebra 0. The root a^ = 5 — 6 is the difference of imaginary 

o 

root 5 and 9 - highest root of the algebra 0. Note that root multiplicity 
mult (a) for affine Lie algebra can be greater than one. 

Subalgebra b+ C g spanned by the generators H\E°' for positive roots 
a G A+ is called Borel subalgebra. 

Parabolic subalgebra p/ D b+ contains Borel subalgebra and is generated 
by some subset of simple roots {aj : j G /, / C {1 . . . r}}. It is spanned by the 
subset of generators {W}U{E°' : a G A+jUJE"" : a G A+, a = Eje/^i"i}- 

Regular subalgebra a C is determined by the root system A^ with the 

6 



set of simple roots {/3j, i = 1, . . . , Va} which is a subset of set of simple roots 
{ai, . . . , ttr-} U {0} with the addition of highest root. 

The Weyl group Wg is generated by reflections {si : f)* — )■ [)*} correspond- 
ing to simple roots {aj}: 

Si- X = \- -Oi (10) 

{ai,ai) 

Root system and characters of representation are invariant with respect to 
the action of Weyl group. Root system can be reconstructed from the set of 
simple roots with the action of Weyl group. 

Weyl groups are finite for finite-dimensional Lie algebras and finitely- 
generated for affine Lie algebras. 

Consider an action of element SaSa+s of Weyl group of afSne Lie algebra 
Q for a simple root of underlying finite-dimensional Lie algebra g. Using the 

definition (10) it is easy to see that SaSa+s-X = '^ + 7^^«+ ( 2f\s) + (•^' '^) ) '^• 
So Weyl group can be presented as semidirect product of Weyl group Wg of 
Q and translations corresponding to the roots of Q. 

Weyl group element can be presented as the product of elementary reflec- 
tions in multiple ways. Number of elementary reflections in shortest sequence 
representing element w G Wg is called length of w and denoted l{w). We also 
use the notation e{w) = (—1 )'(*") for parity of number of Weyl reflections 
generating w. 

Fundamental domain C for the action of Weyl group Wg on P)* is deter- 
mined by the requirement C, E C -x^ {(,,C(i) > for all simple roots aj. It is 
called main Weyl chamber. 

The Cartan matrix A is defined by products of simple roots 

A, = 'P^ (11) 

(«i,«i) 

and can be used for compact description of Lie algebra commutation relations 
in Chevalley basis [18], [22], [23]. 



The form ( 11 ) induces the basis dual to the simple roots basis. It is called 



the fundamental weights basis. We denote its elements by w,: 

{ui, aj) = \ = 5ij (12) 



For finite-dimensional Lie algebra there are r fundamental weights, i = 
1, . . . , r. For affine Lie algebra we have additional fundamental weight uq = 
A, (A, 6) = 1, (A, A) = {6, 6) = 0. Other fundamental weights are equal to 
Ui = a'^Xo+ Ui, where Ui is the fundamental weight of finite-dimensional Lie 
algebra q. 

Weyl vector is equal to the sum of fundamental weights p = ^jCUj. It is 
of importance for the construction of algebra modules. 

2.3. Highest weight modules 

We consider finitely-generated g- modules V such that V = 0fg(,. V^, 
where each V^ is finite-dimensional and there exists finite set of weights 
Ai, . . . As which generates weight system of V, i.e. if dimV^ ^ then C, = 
^i " J2k=i,...,r^kak where Uk G Z+ (See [21], [25]). 

Highest weight module V^ contains one highest weight /i, all other weights 
are obtained by the subtraction of linear combination of simple roots A = 
/i — riiai — ■ ■ ■ — UrCiri nk € Z+. 

The most simple type of highest weight modules is the Verma module 
Af^ who's space can be defined as the space 

M^ = f/(S) ® D^\b+), (13) 

U{b+) 

with respect to the multiplication in [/(g), where means that the action 

U(b+) 

of elements of [/(b+) "falls through" the left part of tensor product onto the 
right part. Here b+ is Borel subalgebra, D'^(b+) is a representation of b+ 
such that D{E'^) = 0, D{H) = fi{H) for any positive root a. Elements of Q 
act from the left and we should commute all the elements of b+ to the right, 
so that they can act on the space Dx{b+). 

Weight multiplicities in Verma module can be found from the Weyl char- 
acter formula 

pM pM 

chMf" = = (14) 

Here we have used the Weyl denominator identity 

R = H [1- e-)'""'*(") = J2 <w)e''^-P, (15) 

aeA+ toSVK 



and e{w) := det (w) is equal to the parity of number of Weyl reflections 
generating w. 

Verma module M^^ has the unique maximal submodule and the unique 
nontrivial simple quotient L^ which is an irreducible highest weight module. 

Irreducible highest weight modules have no non-trivial submodules. Weyl 
character formula for irreducible highest weight modules is 

Thus the character of an irreducible highest weight module can be seen as 
the combination of characters of Verma modules. ( This fact is a consequence 
of the Bernstein-Gelfand-Gelfand resolution flMlEIl, see also 



Construction of generalized Verma modules is analogous to (13), but 
representation of Borel subalgebra is substituted by the representation of 
parabolic subalgebra p/ D b+ generated by some subset {a/} of simple roots 
Jc{l,...,r}: 



Ml' = U (0) ®u(,,) ^p. 



— Q\irLult(cj) 



chMf = — chL^^. (17) 



Denote a formal element Ri := naeA+\A+ (^ ^ ^ °)™'^ • Then the char- 
acter of generalized Verma module can be written as 

1 

Ri 

We can use Weyl character formula to obtain recurrent relation on weight 
multiplicities which can be used for calculations [281 129] . 

For irreducible highest-weight modules recurrent relation has form 

^i = - ^ e(«;)m5_(^(p)_p) + ^ e(w)(5(^(^+p)_p),5. (18) 

w<^W\e w&W 

Formulae for Verma and generalized Verma modules differs only in second 
term on the right-hand side. In case of Verma module it is just (5^^^. For 
generalized Verma module summation in second term on the right-hand side 



(18) is over the Weyl subgroup generated by reflections corresponding to 
roots {a/}. 

Other recurrent formula can be obtained from the study of Casimir ele- 
ments action on irreducible highest weight modules [18] : 

2 

^^ = (n + oV-i\ + o)^ 5Z H^^ + ^"' ")^A+fca. (19) 

9 



It is called the Freudenthal multiplicity formula. Note that it is applicable 
only to irreducible only to irreducible modules. 



We discuss the use of formulae (18) and (19) for computations in section 

m 

Now consider an algebra q and a reductive subalgebra a C g. Simple roots 
Pi of the subalgebra a can be presented as linear combinations of algebra roots 



r„ 



Each 0-module is also an a-module, although L^ is not irreducible as 
a-module. It can be decomposed into the sum of irreducible a-modules: 

L',=®b^.K (20) 

V 

Coefficients in this decomposition are called the branching coefficients. 

It is possible to calculate branching coefficients by the construction and 
repeated subtraction of characters of modules L\. This traditional approach 
has serious limitations especially in case of affine Lie algebras. We discuss 
them in the end of section IH 

Now we describe alternative approach which is based upon recurrent re- 
lation on branching coefficients. But before we proceed to this recurrent 
relation we need several definitions. 

For a subalgebra a C we introduce a subalgebra a^. Consider the root 
subspace \)\^ orthogonal to a, 

and the roots (correspondingly - positive roots) of g orthogonal to a, 

A„^ : ={/3GAg|V/iG[)„;/3(/i)=0}, (21) 

A+^ : ={/3+GA+|V/iG()„;/3+(/i)=0}. 

Let W^^ be the subgroup of W generated by the reflections w^ with the roots 
/3 G A+^ . The subsystem A^^ determines the subalgebra a_L with the Cartan 
subalgebra P)a^. 

Cartan subalgebra can be decomposed in the following way: P) = f)ti © 

We also introduce notation 

aI:=a_L©[)_L a:=a©f)_L. (22) 

10 



For a and a± we consider the corresponding Weyl vectors, Pa and p^j^ and 
compose the so called "defects" Va and "Dnx ^f the injection: 

^a := Pa - T^aP, "^ax := Pax - '^a^P- (23) 

For /i G P^ consider the linked weights {{w{p + p) — p) \w E W} and 
their projections to h^ additionally shifted by the defect — "Pax- 

Pa^iw) ■.= 1Ta^[w{p + p) - p]-'Da^, W E W. 

Among the weights {pai_ (w) \w G W} one can always choose those located 
in the fundamental chamber Ca^ ■ Let U be the set of representatives u for 
the classes W/Wa^ such that 

U:={ueW\ p,A'>^)^CZ_] . (24) 

Thus we can form the subsets: 

p^iu) ■.= -nz[u{p + p)- p]+V^^, ueU, (25) 

and 

Pax («) ■= ^ax HP + P)- p]-'^a^, U E U. (26) 

Notice that the subalgebra a± is regular by definition since it is built on 
a subset of roots of the algebra q. 

Denote by kf' signed branching coefficients. If ^ G Ca is in main Weyl 

chamber kf = bf' otherwise kf = ^iw)bl^,l, \_ where w G Wa is such 

that w{^ + Pa) - Pa^ Wa- 

Now we can use the Weyl character formula to write a recurrent relation 
[30] for signed branching coefficients k/' corresponding to an injection a "^ g: 

The recursion is governed by the set T^^g called the injection fan. The latter 
is defined by the carrier set {'C}a_>n fo^' the coefficient function s{^) 

{^}~^g := {e G P~MO ^ 0} 

appearing in the expansion 

Yl (1 - g-^H") --!*(")— -it»(-^") = _J2 s(7)e-^; (28) 

oeA+\A+ 76Pb 

11 



The weights in (Os-s-n ^^^ ^^ ^e shifted by 70 - the lowest vector in {^} 
and the zero element is to be eliminated: 



a-s>g 



{e-7o|eG{0}\{0}. 



(29) 



The formula (18) is a particular case of recurrent relation for branching co- 
efficients (27) in the case of Cartan subalgebra a = f}g. 



If the root system of a± is generated by some subset of q simple roots 
ai, . . . ,ar then the recurrent relation (27) is connected with the generalized 
Bernstein-Bernstein- Gelf and resolution for parabolic Verma modules [21] ■ 

Another particular case of this formula is connected with tensor product 
decomposition. Consider the tensor product of two irreducible g-modules 
L^ ® L'' . It is also a g-module but not irreducible in general. So 



L'^®L'^ = 0/; 



7 ^ 



(30) 



The coefficients fi^'^ are called fusion coefficients. The problem of computa- 
tion of fusion coefficients is equivalent to branching problem for the diagonal 
subalgebra q C g^Q (see [32]). So our implementation of recurrent algorithm 



can be used to decompose tensor products (See Section 5.1). 

In the case of affine Lie algebras g, a the multiplicities rrii, and the corre- 
sponding branching coefficients b^, can be regarded as the coefficients in power 
series decomposition of string and branching functions correspondingly: 



n=0 j 

00 

bu{q) = ^ K-nsq", ^ = ^ CjUj, Cj > 



(31) 
(32) 



n=0 



String and branching functions have modular and analytic properties which 
are important for conformal field theory, especially in coset models and the 
study of CFT on higher genus surfaces [33], [13], [I2], [3^ . 

AfHne.m calculates weight multiplicities and branching coefficients for 
affine Lie algebras up to some finite grade. We present examples of computa- 



tions in Sections 5.3 , 5.4 Now we proceed to the description of datastructures 



and algorithms implemented in AfRne.m. 
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3. Core datastructures 

Having introduced necessary mathematical objects, problems and rela- 
tions we now describe the related datastructures of AfRne.m. Although 
Mathematica is untyped language it is possible to create structured objects 
and do type checks with patterns [35], [36] . 

3.1. Weights 

Weights are represented by two data-structures: finiteweight for finite- 
dimensional Lie algebras and affineweight for affine. 

Internally the finite weight is a List with the Head f initeweight, its components 
are the coordinates of the weight in the orthogonal basis of Bourbaki [23] . 

Affine weight is an extension of a finite weight by supplying it with level 
and grade coordinates. There is a set of functions defined for finite and affine 
weights. The complete list can be found in online help of the package. The 
most important are the definitions of addition, multiplication by number 
and scalar product (bilinear form) for weights. These definitions allow to use 
traditional notation when working with AfRne.m: 

w=make Finite Weight [{1 ,0 ,3}] ; 
v=niake Finite Weight [{3 ,2 , 1 }] ; 
2*w+v^niakeFinit6Weight [{5 ,2 ,7}] 
w . v==6 

The use of orthogonal basis in the internal structure of weights allows us 
to work with weights without complete specification of root system which is 
useful for study of branching, since roots of the subalgebra can be specified 
by hand. 

3.2. Root systems 

To specify an algebra of finite or affine type it is enough to fix its root 
system. Root systems are represented by two datatypes finiteRootsystem and 
affineRootsystem. The latter is an extension of the former. We offer several 
different constructors for these datastructures. It is possible to specify the 
set of simple roots explicitly, for example to study the subalgebra B2 C B^ 
we can use the definition 

b2b4=makeFiniteRootSystem [ { {1,-1,0,0}, {0,1,0,0} } ] 

There are constructors for root systems of simple finite-dimensional Lie al- 
gebras: 

b2=make S imp leRoot System [B , 2] 
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We use typographic features of Mathematica frontend to offer mathematical 
notation for simple Lie algebras: 

B2 = makeFiniteRootSystem [ { {1, —1}, {0, 1} }] 

Non-twisted affine root systems can be created as afiine extensions of finte 
root systems, e.g. 

b2affine = makeAf f ineExtension [b2 ] 

In notebook interface this can be written simply as B2. 

Semisimple Lie algebras can be created as sums of simple: 

Ai(BAi = f initeRootSystem [2 , 2, {f IniteWeight [2 , {1, 0}], f initeWeight [2 , {0, 1}]}] 

Predicate rootsystemq checks if the object is root systems of finite or afiine 
type. 

List of simple roots is the property of root system so it is accessed as 

rs[simpleRoots]. 

We have implemented several functions to get major properties of root 
systems. Weyl vector is given by the function rho[rs_?rootsystemq]: 

In[l] = rho[b2] 

Out[l] = f initeWeight [2 , {3/2, 1/2}] 

Positive roots can be constructed with positiveRoots[rs_?rootsystemQ]. For afiine 
Lie algebra this and related functions return the list up to fixed grade. This 
grade limit is set through the value of rs[gradeLimit] which is equal to 10 by 
default. List of roots (up to gradeLimit) is returned by roots [rs]. Cartan ma- 
trix and fundamental weights are calculated by the functions cartanMatrix and 
fundamentaiweights Correspondingly. 

It is possible to specify a weight of Lie algebra by its Dynkin label 

weight [b2] [1 ,2 ,3] = makeFiniteWeight [{ 2 , 1}] 

The function dyiikiiiLabeis[rs_?rootsystemQ][wg_?weightQ] retums Dyukiu labels of weight 
wg in the root system rs. 

Elements of Weyl group are specified by the set of indices of reflections, 
so element w = S1S2S1 of Weyl group of algebra B2 is constructed with 
weyiGroupEiement[b2][i,2,i]. Tlicu it cau bc applied to the wcights: 

w = weylGroupElement [b2][l ,2 ,1]; 

w @ makeFiniteWelght [{ 1 ,0}] = makeFiniteWeight [{ - 1 ,0}] 



Computation of lexicographically minimal form [T0l|37j for Weyl group el- 
ements can be conveniently implemented through the use of pattern-matching 
in Mathematica. In [38j rewrite rules for simple finite dimensional and afiine 
Lie algebras are presented as Mathematica patterns. Our presentation of 
Weyl group elements is compatible with the code of [3H] 
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In[l] = <<A3reduce; 

reduce [ s [1 , 2 , 1 ,2 , 1 ,3 , 2 , 1 , 1] ] 
Out [1] = s[2 , 3, 2] 

In[2] = (weylGroupElement [A3] 00 reduce [ s [ 1 , 2 , 1 , 2 , 1 , 3 , 2 , 1 , 1 ] ] ) weight [A3] [- 1 , -2 , - 1] 
Out[2] = f initeWeight [4 , {-2, 2, 1, -1}] 

In[3] = dynkinLabels [yl3] [Out [2]] 
Out [3] = {-4, 1, 2} 

3.3. Formal elements 

We represent formal characters of representations by special structure 
formaiEiement. Tliis structure is a hash-table implemented with DownVaiuos. The 
keys are weights at the exponents and the values are corresponding mul- 
tiplicities. inakeFormalElement[{7i,...,7„},{mi,...,m„}] CrcatCS data-StrUCturC which 

represents the element X^ILi "^i^''^' of formal algebra £. The operations of 
formal algebra £ are implemented for formaiEiement data-type: formal elements 
can be added, multiplied by number or exponent of weight. There is also 
multiplication of formal elements but no division. 

In[l] = makeFormalElement [{ makeFiniteWeight [{ 1 , 1 }] , makeFiniteWeight [{0 ,0}] } , { 1 ,2}] * 

(2 * Exp [makeFiniteWeight [{ 1 ,0}]] * 

makeFormalElement [{ makeFiniteWeight [{1 , 1 }] , makeFiniteWeight [{0 ,0}]} ,{1 ,2}]); 
In [2] = In [1] [weights] 
Out[2] = {f initeWeight [2 , {1, 0}], f initeWeight [2 , {2, 1}], f initeWeight [2 , {3, 2}]} 

In [3] = In [ 1 ][ multiplicities ] 
Out[3] = {8, 8, 2} 

3.4. Modules 

AfRne.m can be used to study different kinds of modules, i.e. Verma 
modules, irreducible modules and parabolic Verma modules . We need 
datastructure module to represent generic module of Lie algebra q. Module 
properties can be deduced from its set of singular weights using Weyl char- 



acter formulae (14), (15), (17), (16). Set of singular weights can have Weyl 
symmetry. It can be symmetry with respect to Weyl group Wq or with re- 
spect to some subalgebra Wa as in the case of parabolic Verma modules. 
Then it is possible to study only main Weyl chamber Ca- To use this sym- 
metry generic constructor for module datastructure accepts several parameters 

makeModule[rs_?rootSystemQ][singWeights_formalElement,subs_?rootSystemQ|emptyRootSystem[],limit:10. 

Here rs is root system of Lie algebra 0,singweights is the set of singular weights, 
subs is root system corresponding to Weyl group Wa which is the symmetry of 
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the set of singular weights. Parameter limit hniits computation for infinite- 
dimensional modules such as Verma or parabolic Verma. There are several 
specialized constructors for different types of highest weight modules: 

vin=makeVermaModule [B2 ] [{ 2 , 1 } ] ; 

piii=makeParabolicVermaModule [-B2 ] [ weight [B2][2 ,1] ,{1}]; 
im=makeIrreducibleModule [B2 ] [ 2 , 1 ] ; 
GraphicsRow[textPlot/@{ini,vm,pm}] 
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As we already stated properties of the module are encoded by its singular 
element. Function singuiarEiemeiit[m_moduie] returns singular element of module 
as formaiEiement datastructure. Character (up to limit for (parabolic) Verma 
modules) is returned by function character [m.moduie]. Direct sum of modules is 
module and we use natural notation 

inil=niakeIrreducibleModule [B2 ] [weight [B2 ] [2 , 1 ] ] ; 
im2=niakeIrreducibleModule [B2 ] [weight [B2 ] [1 , 2 ] ] ; 
textPlot f iml© im2l 
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Tensor product is also implemented but only for finite-dimensional Lie 
algebras, since tensor product of afiine Lie algebra modules leads to rich new 
structures f39l HOl E] which are out of the scope of present paper. 
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textPlot [makelrreducibleModule [Ai ] [5] Cg) makelrreducibleModule [Ai ] [3] 
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4. Computational algorithms 



As we have already stated in section 2.3 there exist two recurrent re 



lations which can be used to calculate weight multiplicities in irreducible 
modules. Both algorithms proceed in the following way to calculate weight 
multiplicities: 

1. Create the list of weights in main Weyl chamber by subtraction of 
all possible combinations of simple roots from the highest weight (e.g. 
for finite-dimensional algebra subtract ai from fi while inside C, then 
subtract 02 from all the weights already obtained etc). 

2. Sort the list of weights by their product with Weyl vector. 

3. Use a recurrent formula. If the weight required for recurrent computa- 
tion is outside the main chamber use Weyl symmetry. 

The difference in performance of algorithms is in the number of previous val- 
ues required to get the multiplicity of weight under consideration. For Weyl 



formula-based recurrent relation ( 18 ) it is constant and equal to the number 



of elements in Weyl group (if we are far from the boundary of representation 
diagram). When Freudenthal formula (19) is used number of previous val- 



ues grows with the distance from the external border of representation. So 
Freudenthal formula is faster if the weight is close to the border or the rank 
of the algebra and the size of Weyl group is large [8j . Note that Freudenthal 
formula is valid for irreducible modules only, so it can not be used to study 
(generalized) Verma modules. 

We have made some experiments with our implementations of Freuden- 
thal formula and formula ( 19 ) and have got Figure [II which depicts depen- 
dence of computation time on number of weights in module. 
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Figure 1: Running time of algorithms based on Freudenthal formula (19) (dashed) and 



recurrent relation ( 18 1 (solid) with the number of weights in C for calculation of multi- 



plicities in representations of -62- 



In the calculation of branching coefficients use of Freudenthal formula 
requires full construction of formal characters of algebra representation and 
all the representations of subalgebra. It is impractical when the rank of 
algebra and subalgebra are big, for example for maximal subalgebra. 

Alternative algorithm which was presented in the paper [20! proceeds as 
follows: It contains the following steps: 

1. Construct the root system A^ for the embedding a — ?■ g. 

2. Select all positive roots a G A+ orthogonal to a, i.e. form the set A+ . 

3. Construct the set T^^^- Relation ^ defines the sign function 5(7) and 
the set $ocB where the lowest weight 70 is to be subtracted to get the 
fan: T^^g = ^ - 7o|e e *acj \ {0}. 



4. Construct the set \I'('') 
for the 0-module L^^^. 

5. Select the weights {/ij^ 



{w(yU + p) — p; u; G W} of singular weights 



\w 



n 



a± 



\w 



(p + p) - p] - Va^ eC^}. Since 



the set A+^ is fixed we can easily check wether the weight fi^ (w) be 
longs to the main Weyl chamber C^ (by computing its scalar product 
with the fundamental weights of a']_). 
6. For the weights pj-^ (w) calculate dimensions of the corresponding mod- 



ules, dim ( L 



'H 



the singular element \E'/|^ 



using the Weyl dimension formula and construct 

(m) 



(o,a±)- 



7. Calculate the anomalous branching coefficients using the recurrent re- 



lation (27) and select among them those corresponding to the weights 



in the main Weyl chamber Ca- 

We can speed up the algorithm by one-time computation of the represen- 
tatives of the conjugate classes W/W^^ ■ 

Consider the regular embedding B2 C -B4. In this case fan consists of 24 
elements. In order to decompose B^ module we need to construct the subset 
of singular weights of the module which projects to the main Weyl chamber of 
subalgebra -82- Full set of singular weights consists of 384 elements. Required 
subset contains at most 48 elements. So time of the construction of required 
subset is negligible if the number of branching coefficients is greater than 
that. We may estimate the total number of required operations for the 
computation of branching coefficients as the product of number of elements 
in main Weyl chamber of subalgebra with non-zero branching coefficients 
and number of elements in fan. In the case of direct algorithm we need 
to compute the multiplicities for each module in the decomposition, so the 
number of operations grows faster than square of number of elements in main 
Weyl chamber of subalgebra with non-zero branching coefficients. 

To further illustrate this performance issue we include the Figure [2] where 
we show the time required for computations of branching coefficients for 
53 C S4. 

5. Examples 

In this section we present some examples of computations available with 
AfRne.m with the code required to produce these results. 

5.1. Tensor product decompositon for finite-dimensional Lie algebras 

Computation of fusion coefficients for the decomposition of tensor prod- 
uct of highest-weight modules to the direct sum of irreducible modules has 
numerous applications in physics. For example, we can consider spin of 
composite system such as atom. Another interesting example is integrable 
spin chain consisting of N particles with the spins living in some representa- 
tion L of Lie algebra q with g-invariant Hamiltonian H, describing nearest- 
neighbour spin-spin interaction. In order to solve such system, i.e. find 
eigenstates of Hamiltonian, we need to decompose L*^^ into the direct sum 
of irreducible g-modules of lower dimension and diagonalize the Hamiltonian 
on these modules. 

19 



14 



12 



10 



10 



15 



20 



25 



30 



35 



Figure 2: Running time of algorithms based on Freudenthal formula (19) (dashed) and 



recurrent relation ( 27 ) (solid) with the number of weights in C for calculation of branching 
coefficients for B^ C B4. 



For fundamental representations of simple Lie algebras it is sometimes 
possible to get analytic result for the dependence of decomposition coeffi- 
cients on A^ (See [32] )• Our code give numerical values and can be used to 
check this analytic results. 

Consider tensor power of B2 first fundamental representations (L^^'^l) 
Decomposition coefficients are just branching coefficients for tensor product 
module to the diagonal subalgebra B2 <Z B2 ® B2 ® B2 ® B2. So following 
code calculates these coefficients: 

fm = makelrreducibleModule [,82 ] [ 1 1 0] ; 
tp = ((fm® ]fm)(g) fm)®]fm; 
subs = makeFiniteRootSystem [ 

{1/4*{1, -1, 1, -1, 1, -1, 1, -1}, 
l/4*{0, 1, 0, 1, 0, 1, 0, 1}}]; 
be = branching [tp , subs]; 
{bc[#], dynkinLabels [ subs ] [#] } & /a bc[weights] 

It produces list of highest weights and tensor product decomposition coeffi- 
cients: 

{{1, {4, 0}}, {3, {2, 2}}, {0, {3, 0}}, 
{2, {0, 4}}, {3, {1, 2}}, {6, {2, 0}}, 
{6, {0, 2}}, {1, {1, 0}}, {3, {0, 0}}}] 

Returning to the problem of spin chain Hamiltonian diagonalization we 
can see that instead of diagonalizing operator in space of dimension 625 we 
can diagonalize operators in spaces of dimensions 55, 81, 30, 35, 35, 14, 10, 5, 1. 
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5.2. Branching and parabolic Verma modules 

We illustrate generalized BGG-resolution with the diagrams of G2 parabolic 
Verma modules which appear in the decomposition of irreducible module 



ch (L^) = J2 e^»("^e(M)chM; 



f^a I (") 



(33) 



ueu 



Character of L^^'^J is presented in Figure [3| characters of generalized Verma 
modules in decomposition (33) are shown in Figure |4J Characters in upper 
row appear in (33) with positive sign and in lower row with negative. 
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Figure 3: Character of irreducible G2-niodule L^^'^l 

5.3. String functions of affine Lie algebras and CFT models 

String functions can be used to present the formal character of affine Lie 
algebra highest weight representation. They have interesting analytic and 
modular properties [HI [321 1^ - 

AfRne.m produces power series decomposition for string functions. Con- 
sider affine Lie algebra s/(3) = A2 highest weight module L^^'^'^\ To get 
string functions we can use the code: 
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Figure 4: Character of G2 generalized Vernia modules appearing in decomposition of 
-L'^'^1. Parabolic Verma modules in upper row appear in decomposition with positive sign, 
in lower row - with negative. 



str ingFunctions [A2 ,{1 ,1 ,2}] 
{{0, 0, 4}, 

2q + 10g2 + iOq'i + 133g4 + 393^5 _|_ io84g*' + 2760g^ + 6632g» + 15214g9 + 33508gl'' } , 
{{0, 3, 1}, 

2q + 12g2 + 49gS 4, iQQg4, ^ 494^5 _)_ 1340^6 ^ 3387^7 + gosGcj* + 18415<?8 + 40302gl0 } , 

{{1, 1, 2}, 

l + 6q + 27<j2 + 96qi + 298<?'' + 8369"' + 21739" + 5310g'^ + 12341g** + 27486g9 + 59029(j"' } , 
{{2, 2, 0}, 

1 + 8<2 + 35<j2 + I24g3 + 379g-* + 1052g5 + 2700g*' + 6536g7 + 15047(?« + 33248g9 + 70877g"' } , 
{{3, 0, 1}, 

2 + 12g + 49g2 + 166<?''' + 494g'' + 1340g5 + 3387?'' + 8086g^ + 18415g* + 40302g9 + 85226<?"'} 

Similarly for affine Lie algebra G2 we get 

str ingFunctions [G2 ,{1 ,1 ,0}] 
{{2, 0, 0}, 

1 + 813 + 37(j2 + 138(j» + 431q* + 1227g5 -f 3208(3*^ + 7901g'^ } , 
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{{0, 0, 1}, 

3g + 18g^ + 735^ + 2'i7q^ + 736q^ + 2000g« + 5070g'' } , 
{{1, 1, 0}, 

1 + 7q + 32g2 + Il7g3 _,_ 379^4 ^ I055g5 + 2780(j» + GSSOq"^ } , 
{{0, 2, 0}, 

3g + 15(j2 + 63?=* + 210?"' + 633g5 + 1725g^ + 4407g''} 

5.4- Branching functions and coset models of conformal field theory 

More complex models of CFT can be obtained from cosets G/A corre- 
sponding to the embedding a C 0. These models can be studied as gauge 
theories ^E^. 

Branching functions of the embedding a C g are the partition functions 
of CFT on the torus (see [I3]). 

As first example we show computation of branching functions for the 
embedding Ai — )■ B2 up to tenth grade: 

branchingFunct ions [B2 ,makeAffineExtension[inakeFiniteRootSystein[{{l, 1}}]], {1, 1, 1}] 

{{3, 0}, 

2 + 14ij + 52g2 + 154(j3 + 410?"' + 994i}5 + 2248ij^ + 4832ij'' + 9934ij** + 19680g8 + 37802iji'' } , 
{{2, 1}, 

4 + 20g + 72q2 + 220g3 + 584g4 + 1424g5 + 3248^'' + 7012g7 + 14488g* + 28844(j9 + 55616qi0 } , 
{{0, 3}, 

4g + 20^2 + 68?=* + 200<j'' + 516g5 + 1224g6 + 2736g'' + 5808g* + 11820g8 + 23236g"^ } , 
{{1, 2}, 

2 + 14g + 54g2 + 168(j3 + 462g"' + 1148^5 + 2656?'' + 5812g^ + 12130(3* + 24358<79 + 47328(?10 } 

Another example is computation of branching functions for the regular 
embedding i?2 C Cs: 

sub=makeAf f ineExtension [parabolicSubalgebra [C3] [2 ,3]] ; 
branchingFunct ions [Cs , sub , {2, 0, 0, 0}] 

{{0, 1, 0}, 

2q - 20g3 + 24.q'^ + 82g5 _ 320q6 + lOSq"^ } , 
{{1, 0, 0}, 

1 - g - 8g2 + 19g3 + Wq^ - 156g5 + 205(j'^ + GiOq"^ } , 
{{0, 0, 1}, 

q-5q'-^ + 7q^ } 



6. Conclusion 

We have presented the package AfRne.m for computations in represen- 
tation theory of finite-dimensional and afiine Lie algebras. It can be used to 
study Weyl groups, roots systems, irreducible, Verma and parabolic Verma 
modules of finite-dimensional and affine Lie algebras. In present paper we 
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have also discussed main ideas used for implementation of the package and 
described most important notions of representation theory required to use 
AfRne.m. 

We have demonstrated that recurrent approach based upon Weyl charac- 
ter formula is not only useful for calculations but also allows to see connection 
with (generalized) Bernstein-Bernstein- Gelf and resolution. 

Also we have presented examples of computations with this package con- 
nected with problems of physics and mathematics. 

In future versions of our software we are going to treat twisted affine Lie 
algebras, extended affine Lie algebras and provide more direct support for 
tensor product decomposition. 
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Appendix A. Software package 



The package can be freely downloaded from Ihttp : //github . com/naa/ 
|Aff ine( To get the development code use the command 

git clone git :// github . com/naa/ Aff ine . git 

Contents of the package: 

Affine/ 
demo/ 
demo . nb 
paper . nb 
doc/ 

figures/ 
timing.pdf 

branching- timing . pdf 
irrep-sum.pdf 
irrep-verma-pverma . pdf 
G2-irrep.pdf 



root folder 

demonstrations 
demo notebook 
code for the paper 
documentation folder 
figures in paper 

diagram showing performance 

for branching coefficients 
sum of B2 irreps 
irrep, Verma, (p)Verma for B2 
irrep for G2 
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G2-pverma.pdf 
tensor-product .pdf 

bibliography . bib 

paper.pdf 

paper.tex 

TODO . org 
src/ 

aff ine.m 
tests/ 

tests. m 
README . markdown 



parabolic Verma for G2 
tensor product of Al-modules 

bibliographic database 

present paper 

paper source 

list of issues 
source folder 

main software package 
unit tests folder 

unit tests 
installation and usage notes 
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